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Abstract 

We study the generahzed harmonic oscillator which has both the position-dependent mass and 
the potential depending on the form of mass function in a more general framework. The explicit 
expressions of the eigenvalue and eigenfunction for such system are given, they have the same 
forms as those for the usual harmonic oscillator with constant mass. The coherent state and the 
its properties for the system with PDM are also discussed. We give the corresponding effective 
potentials for several mass functions, the systems with such potentials are isospectral to the usual 
harmonic oscillator. 
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I. INTRODUCTION 



Recently, the study of quantum systems with position-dependent mass (PPM) has at 
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321 ] ■ The models have played important 
roles in description of the electronic properties of semiconductors [1], quantum dots {4], 
liquid crystals 5| and so on. In the theoretical works on the system with PDM, the main 
concern was the exact solutions of the corresponding Schrodinger equation or its relativis- 



tic generalizations such ^_as K 
transformation method 
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ein-Gordon equation and Dirac equation. Using coordinate 



35], supersymmetric quantum mechanics [Sg, l37| and other 



methods, many solvable potentials of the wave equations for different mass functions have 
been obtained. These potentials can be classified according to the forms of eigenfunctions 
of the wave equations, which show some similarities between systems with constant mass 



and with PDM 



3^- 



It is well known that the harmonic oscillator with constant mass is a very important 
system in quantum mechanics 



38!]. Apart from its wide applications in problems of con- 
densed matter, atomic, nuclear and elementary particle physics, the harmonic oscillator also 
has connections with many interesting methods analytically solving Schrodinger equation 
such as factorization method, supersymmetric quantum mechanics, algebraic method. So, 
it is natural to generalize the usual harmonic oscillator and related studies to the case of 
position-dependent mass. There exist some works in this aspect by using point canoni- 



cal transformation and Lie algebraic approach[25|, |26|, |27|. However, these investigations 
only involve the specific potentials and mass functions. In this paper, we consider the so- 
called generalized harmonic oscillator(GIIO) based on the following two requirements: the 
Hamiltonian for the system with PDM can be factorized in terms of the rising operator and 
lowering operator which are generalization of the creation and destruction operators for the 
usual harmonic oscillator, respectively; the rising operator and lowering operator satisfy the 
same form of commutation relation as that for the creation and destruction operators. With 
these requirements, we put the discussion of GHO in a more general framework. We will 
see that many systems with PDM rather than with specific mass functions have the same 
energy spectrum and same form of eigenfunctions as those of the usual harmonic oscillator. 
The paper is organized as follows. In section II, the condition for the Hamiltonian with 



2 



PDM to be that for GHO is derived, the eigenproblem of GHO is solved using the operator 
method. Some studies such as the coherent states related to the harmonic oscillator are 
generalized to the system with PDM and their properties are discussed. In section III, for 
some mass functions, the corresponding potentials are given. In the last section, we will 
make some remarks. 



II. GENERALIZED HARMONIC OSCILLATOR AND ITS SOLUTION 

When the mass of a particle depends on its position, the mass and momentum operators 
no longer commute j^, so there are several ways to define the kinetic energy operator. We 
start with the kinetic energy introduced by von Roos[2], the Hamiltonian with position- 
dependent mass M(r) and potential energy V{r) reads 
1 



+ V^(r) 



(1) 



H = ^ ^M{r)''pM{r)"pM{ry + M(r)^pM(f)''pM(f)" + V{r) 

where mo is a constant mass and m(r) is a dimensionless position-dependent mass M(r) = 
mom{r), a, f3, 7 are parameters and satisfy the condition ofa + /3 + 7 = — 1. Using natural 
units (mo = h = 1) and only considering one- dimensional system, we can rewrite Hamilto- 
nian ([1]) as 

1 d 1 d 1 d"^ m' d , . 



where the effective potential is 



= Vix) + + ^[«(« + P + + 1]^, (3) 



andm'(x) = ^,m"(a:)- 

There are many debates for the choice of the parameters a, /5, 7. Morrow and 
Brownsteinj^ have shown that = 7 based on the comparison between the experimen- 
tal results and the analytical solutions of some models. In the following, we use this kind 
of condition for the parameters as a basic fact. We will introduce the lowering and rising 
operators A and so that the Hamiltonian ([2]) can be factorized in terms of these operators 
and its solution be obtained by using the operator method. For this purpose, we define the 
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lowering and rising operators A and as follows 17|, respectively 



V2 



V2 



ax 



ax 



where the function will be determined by the requirement that 

[i,i+] = l. 



(4) 
(5) 

(6) 



The constrain ([6]) makes our discussion be restricted to GHO, so u{x) in Eqs.(jl]) and ([5]) is 



not the so-called superpotential, which is different from that in 
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25 



26| 



Using Eqs.(jl]) and ([5]) and performing some calculations, we have 



^^'{x) 



'mix] 



4m(x 



m"{x) 3 / m'(x) 



m{x) 2 \ m{x) 



(7) 



It follows from Eqs.(l6]) and ([7]) that fi{x) and (3 are given by the following relations, respec- 
tively 

(3 = --, fi{x) = / ^/m{x)dx. (8) 
Now, Eqs.(jl]) and ([5]) can be rewritten as follows. 



V2 



A' 



1 
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m(x) ^—m(x) 4 
ax 

, , 1 d 



(9) 



It is easy to see that the operators A and will become the destruction and creation 
operators for the harmonic oscillator respectively when the mass m{x) is a constant. 

If we assume that the potential function V{x) in Eq.(l2]) depends on mass m{x) and is 
given by 

(10) 



V{x) = -Ux)]\ 



then the Hamiltonian (12]) can be rewritten as 



H = A+A+- = iV+ -, 
2 2' 



(11) 



where N = A'^A may be called number-like operator because it reduces to the number 
operator for the harmonic oscillator when the mass m{x) is independent of x. It is also 



evident that with m{x) = 1 the potential (fTOj) and the Hamiltonian (fTTj) become the potential 
and Hamiltonian for the harmonic oscillator, respectively. In this sense, we call the above 
system with PDM a generalized harmonic oscillator. 

With Eq.([n]), we can easily get the commutators between operators A, ^4+ and N 

[N,A] = -A, [iV,i+] = i+. (12) 



The following procedure solving the eigenproblem of Eq. ffTT]) is similar to that of the har- 
monic oscillator [39]. If we assume that i/Jn{x) is an eigenstate of N with eigenvalue n, that 
is 

Nipn{x) = nipnix), (13) 
then using Eg. (1121) . we have the relation 

nA+Mx) = {n + l)i+^„(x). (14) 

Eq. ffT^ indicates that A'^tpni^) is also an eigenstate of with eigenvalue [n + 1), so A'^ 
can be regarded as a raising operator. By a similar reasoning, A is a lowering operator 
which means that Aipnix) is an eigenstate of N with eigenvalue (n — 1). If we assume the 
eigenstate of has a lower bound, that is there exists a state i/jo{x) satisfying the relation 



m(x) 4 — (m(a;) "^tpolx)) + n(x)tpQ(x) 
ax 



0, (15) 



then all eigenstates of N can be obtained by successive application of the operator A'^ on the 
state ipo{x). The existence of the lower bound is related to the requirement that |A'0„(a;)|^ 
is real and positive. Solving Eq.(|T5l). we get the analytical expression of ipoix) 

i/joix) = [m(x)]3e-^['^(^)]'. (16) 

Now, acting ipo{x) on the left with the operator we have 

/2 

iPi{x) = A+i/joix) = V2fi{x)M^) = ^[m(x)]^e-5['^(")l'i/i(^(x)). (17) 
Similarly , we get 

Mx) = (i+)¥o(x) = {2[^(x)]2 - l}Mx) = ^[m{x)]-*e--2^'^^^^^'HMx)). (18) 
By induction, one can show that all eigenstates of H are given by 

^„(x) =A/;Hx)]ie'^['^(^)l'i/„(Mx)), (n = 0,l,2,...), (19) 



where Mn is the normahzation coefficient, /i(a;) is defined by Eq.([8]), and Hn{^{x)) is the 
Hermite polynomiaL The normahzation of -i/^n means that 

/oo 

°° ('20'! 



rmn 



where //mm = fJ-i—oo) and /imaa; = yu(oo). If /imin = or — cxD and fimax = oo, we can get the 
exphcit expression of J\fn as follows 

1 



I ■ — p=i il-I'min 00) 

^ f (21) 



The requirement for the parameters fimim fJ-max will input constraint on the mass functions. 
However, if we require that ipn{x) is orthogonal for different n, then there is only one choice 

of l-^min 

Using Eqs. (fTT]) and (IT^ . we have 

HM^)= + (n = 0,l,2,...), (22) 

which means that ipn{,x) is the eigenstate of H with eigenvalue En = n + \. In a sense, GHO 
is in the same class as the harmonic oscillator considering that they have the same energy 
spectrum but with different potentials. 

Now, we can calculate the matrix elements {n\A^\n'), {n\A\n') which are the same as 
those for the harmonic oscillator 

1*00 /"CXD 

(n|A+|n') = / ^nA'^'^n' dx = Vn' + 1 / Vn^n'+i dx = Vn/~TTSn,n'+l, 

J-oo ^23) 



— oo 

CO 



{n\A\n') = ipnMn' dx = ^/ n' I i^^i^n'-idx = Vn'6ny-i 



Because the operators A and ^4+ have the same kind of properties of the creation and 
destruction operators, many studies on the latter may be generalized to the former. For 
example, the canonical coherent state is defined to be the eigenstate of the destruction 
operator jiol. If we generalize this definition of the coherent state to the operator A, that is, 

A\z) = z\z), (24) 



6 



where z is a complex number, and similarly we expand \z) in term of the state \z) 



c„|n), then we have 

n=0 



n=0 



(25) 



where the factor e~~ comes from the convention of {z\z) = 1. The completeness relations, 
non-orthogonality for different eigenstates z and z' , and other properties for the coherent 
state of system with constant mass are also hold for the system with PDM. However, due 
to the fact that we cannot express x in terms of the combinations of the operators A and 
as in the case of harmonic oscillator, it is difficult to calculate the matrix elements such 
as {z\x\z), {z\x^\z). From Eqs.([9]), we have 

(26) 



/i(x) 
7r(a;) 



—im 4— m 4 
ax 



where ir is the so-called deformed momentum 
expressions of the matrix elements 

1 



2l| . With Eq. (124|) , it is easy to get the following 



{z\fi{x)\z) 



V2 



[z+z*), {z\f,{xr\z) = -[{z+zr+i], 



{z\t({x)\z) = - — {Z- Z*), {z\Tl{xf\z) 



1], 



(27) 



and 



A/i(x) = ^/ {z\^i{xY\z) - {{z\^i{x)\z)) 
A7r(x) = a/ {z\tc{xY\z) — {{z\t[{x)\z)) 



1 



(28) 



Noticing that ii{x) and 7r(z) are functions of x except that 7r(z) is related to the momentum 
p, Eqs. (1281) show that there exist more quantities with minimum uncertainty in the coherent 
states for systems with PDM than for those with constant mass. 

It should be emphasized that above discussion is the generalization of the method used for 
the usual harmonic oscillator. In this setting, the analytical expressions of the eigenvalues 
and eigenfunctions for GHO are obtained through the so-called Heisenberg algebra generated 
by A^^A^N and the unit operator 1 rather than the Lie algebra sm(1,1) or the nonlinear 
algebra studied in 25|, l26|]. It is also noted that the condition (6) was used to determine the 



superpotential in the framework of supersymmetry quantum mechanics in 25|, |26|] . However, 
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the condition (6) has several functions in our discussions: (i) determining the parameters 
in the Hamiltonian ([2]); (ii) indicating that the potential in the Hamiltonian ([2]) is related 
to the mass function when the Hamiltonian is required to be factorized with the operators 
obeying the condition (6); (iii) showing that the algebraic methods and some related studies 
to the usual harmonic oscillator may be applicable to the systems with PDM. As to (iii), the 
coherent state for the system of PDM discussed above is such an example and it gives a more 
general result on minimum uncertainty than that for the case of constant mass. In addition, 
we also study the orthonormality of the wavefunctions and its possible restriction to the 
form of the mass functions, which is illustrated through examples in the next section. To 
our knowledge, all those mentioned above are not appeared in the literature for the systems 
with PDM. 



III. MASS FUNCTIONS AND POTENTIALS 



In the last section, we have shown that the generalized harmonic oscillators have the same 
energy spectrums and the same form eigenstates as those of the harmonic oscillators for a 
lot of mass functions. In this section, we will give the corresponding potentials for several 
mass functions that satisfy the requirements for the orthonormality of i^nix)- 

Example 1. We first consider the mass function m{x) to be of the form 



mix) 



a + X 



(29) 



1 + X2 

where a > 0. This effective mass was used by many authors in their studies on solvability of 



14 



15 



17 



25 



26|. It is noted that we have 



the systems with PDM 
which corresponds to the harmonic oscillator. 
Substituting Eq. (l29!l into Eq.([8]), we have 

IJ,{x) = X + {a — 1) arctanx, 

and the related potential function given by 

2 , (a - 1) [3x^ + (4 - 2a)x^ 



mix) 



1 for a 



Veffix) = lHx)r+^—^ 



(30) 



(31) 



(a + a;2)4 

which is the effective-mass analogue of the singular oscillator potential. From Eqs. fllQI) and 
f l29|) . we get the eigenstate as follows 



M^) = A/;y^^e-^I^(^)l'i/„(/i(x)), [n = 0, 1, 2, ...). 



(32) 



These results are the same as those obtained b 
SUSY quantum mechanics method and so onjo], [l4 
Examj) 



;he coordinates transformation method, 



15|, 



25|. 



the form 
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e 2. The mass for a particle is exponentially increasing or decreasing which has 



mix) 



(33) 



This kind of mass function may be used in the study of confine energy states for carriers in 
semiconductor quantum well structures I2I, 3- For the mass function (1551) . Eq.([H]) gives 



2 

-62 

a 



(34) 



(a 7^0), 
X, (a = 0), 

and a = is related to m(x) = 1, that is to the case of constant mass. From Eq. fllOl) . we 
have the potential 

' ^ (a^O), 



V{x) 



-X 

2 



[a = 0), 



(35) 



which indicates it has a similar behavior as the mass (133!) for a 7^ 0[ll^. The effective 
potential corresponding to the mass fl5^ can be written as 



Veff{x) = \[f.{x)r~^e-^\ 



The corresponding wavefunctions are 

^„(x) =Ar„et-e-^[^(^')]'//„(Mx)), (n = 0,1,2,...). 



(36) 



(37) 



It is noted that the orthogonality of the eigenstate requires a < for x < and a > for 
x > 0. 

Example 3. We consider another mass which reads as 



m{x) = 1 + tanh(aa;). 



(38) 



where a is a real parameter. The solution of Schrodinger equation with mass fl38l) and some 
smooth potential was studied in Refs.[3] and [8|] by solving Schrodinger equation directly. 
In this case, we get the expression for fi{x) in Eq. ([8]) 



fi{x) 



— ln(e"^ + VlT^), (a ^ 0) , 
a 

x, (a = 0). 



(39) 



The effective potential corresponding to the mass (l38i) is 
and the wavefunctions are 



K//(a:) = -[Mx)f-- ^,,, " (40) 



^^(2:) = Ar„[l + tanh(ax)]3e-5[MW] i/^(^(a;)), (n = 0, 1, 2, ...)• (41) 
Example 4. The mass is a power of x, that is 

m(x) = x", (42) 
where a is a real parameter. For this mass, fi{x) in Eq. ([8]) takes the form 

2 (a ^-2), 



lj{x) = <( a + 2 (43) 
Inx, (a = —2). 

The effective potential and the wavefunctions related to the mass (142|) are, respectively 

Veff = Ih^)? - ^a{3a + A)x-^^^'\ (44) 
^Pn{x) = A/;xte-^IM-)l^if„(;,(:,)), (n = 0, 1, 2, ...)• (45) 
Example 5. The mass is a deformed hyperbolic function of x 

m{x) = sech^(ax), (46) 

where a is a positive parameter. Now, //(x) in Eq. ([8]) is 

2 

/i(2;) = - arctane""". (47) 
The corresponding effective potential is 

Veff = ^H^)? - ^[3cosh(2ax) + 1]. (48) 

The wavefunctions are 

M^) = K[sech{ax)]h--^^>'^^^^" HMx)), (n = 0, 1, 2, ...). (49) 



Example 6. The mass function has the form 27 1 



,2 



= (50) 
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where a is a real parameter and q > 0. When g = 0, (H6!) reduces to that of (H2!) . In this 
case, fi{x) in Eq. ([8]) is 

ll(x) = arctan ( | . (51) 



The effective potential is 

Veff = — arctan 



X 

71 



q X 



(52) 



(53) 



2q 

The wavef unctions are 

^n(x) = AT /^Ie-i[^(^)l'if„(//(x)), (n = 0, 1, 2, ...)• 
\l q + x'^ 

It is noted that mass (l50l) has no analogue of the constant mass for whatever a and which 
is different from the mass functions in the other five examples. Also, the orthogonality of 
the wavefunction for GHO requires that /i(x) have the property //mm = — c>o and ^imax = oo, 
but /i(x) in Eqs. Wf} and fl5T]) do not so. In this sense, the mass functions fj46l) and fl50|) 
should be excluded from a possible choice of the mass, and the wavefunctions in (H9l) and 
fl53|) are only formal. Note that the wavefunctions in ^] are unnormalized. 



IV. REMARKS AND DISCUSSIONS 



In the above sections, we have discussed a special system with PDM, i.e. the generalized 
harmonic oscillator that can be solved using operator method. We obtain the analytical 
expressions of its energy spectrums and eigenstates. For several mass functions, the corre- 
sponding effective potentials are given. In our discussions, we restrict the lowering and rising 
operators to obey the commutation relation ([6]) which is same as that for the creation and 
destruction operators of the harmonic oscillator. With this kind of construction, we have 
restrict the choice of the parameters a,/3,7 in the Hamiltonian ([2]). If we further require 
that the Hamiltonian ([2]) can be factorized into the product of the lowering and rising oper- 
ators, then the potential V{x) is determined by the mass function, which is consistent with 
that from the SUSY quantum mechanics method. Different from the SUSY method, the 
key of our discussions lies in that we can generalize a lot of investigations of the harmonic 
oscillator to GHO. For example, the coherent state in sect ion [IT] for the system with PDM is 
studied with this kind of consideration, and it indeed shows some differences from that for 
the system of constant mass. Also, the orthonormality of the wavefunctions for GHO puts 
the restriction on the choice of the mass functions. 
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